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Abstract 

We investigate the correlators of TrA^A y in matrix models on homogeneous spaces: 
S 2 and S 2 x S 2 . Their expectation value is a good order parameter to measure the 
geometry of the space on which non-commutative gauge theory is realized. They also 
serve as the Wilson lines which carry the minimum momentum. We develop an efficient 
procedure to calculate them through 1PI diagrams. We determine the large iV scaling 
behavior of the correlators. The order parameter shows that fuzzy S 2 x S 2 acquires 
a 4 dimensional fractal structure in contrast to fuzzy S 2 . We also find that the two 
point functions exhibit logarithmic scaling violations. 
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1 Introduction 



One of the most important problems in string theory is to understand the dynamics of D- 
branes. The understanding of such a problem is likely to lead to the ultimate formulation 
of string theory. For example, old matrix models for non critical strings are reinterpreted as 
matrix models of D-branes[l][2]. 

Such a development reinforces our interest in matrix models for critical strings [3] [4]. In 
matrix models, non- commutative solutions are identified as D-branes. On such a solution, 
non- commutative (NC) gauge theory is realized [5] [6] [7]. In string theory, NC gauge theory 
arises on D-branes with B^ v field[8]. More recently the non-commutativity parameter is 
identified with g s a' where g s is the string coupling constant and a' is the inverse string 
tension[9][10]. We therefore suspect that non-commutativity plays a fundamental role in 
string theory [11]. 

In a non-perturbative investigation, it is desirable to work with matrices of finite N. In 
such a setting, branes are inevitably compact. The simplest such objects in matrix models 
are fuzzy sphere and its higher dimensional analogues[12][13][14]. In order to obtain such a 
classical solution in matrix models with finite N, we need to introduce Myers term and its 
generalizations[15][16][17]. It serves as a nonperturbative formulation of NC gauge theories 
on homogeneous spaces G/H. 

We have investigated quantum corrections of NC gauge theory on G/if [18]~[21]. We 
have identified the 't Hooft couplings and large N scaling behavior of these gauge theories 
by power counting arguments. Recently a non-perturbative study has been performed in a 
bosonic matrix model[22]. The large N scaling behavior of the observables are found to be 
in good agreement with perturbative predictions. 

From string theory point of view, the Myers term corresponds to the constant H^ vp 
field. With finite N, most of the solutions are metastable since the single S 2 solution always 
minimizes the action at the classical level. Nevertheless we expect that they are stabilized 
in the large N limit due to the suppression of the tunneling effect. We have indeed found 
that the H^p field vanishes in our large N scaling analysis. We therefore argue that various 
G/H solutions approach non-commutative R d with various dimensionality d in the large N 
limit realizing the formal classical solutions of IIB matrix model. Our investigations concern 
quantum fluctuations around a particular extremum of the effective action and the exploring 
the entire landscape is beyond the scope of this paper. 



1 



Another possibility is to contemplate that such a solution extremizes the quantum ef- 
fective action of IIB matrix model. In fact we have found that a 4 dimensional solution of 
S 2 x S 2 type extremizes the effective action at the two loop level[19][21]. Since its effective 
action is O(N) which is the same order with the space-time volume, we may interpret this 
solution as a D3 brane with a finite tension. It is because the true minimum of the effective 
action is argued to be 0(1) in IIB matrix model [23] [24]. 

In the matrix model construction of NC gauge theory space-time and gauge field are 
unified into the identical matrix degrees of freedom which are identified as the coordinates 
semiclassically. In this sense the gauge invariant operator TrA^A u is a good order parameter 
to measure the shape and extension of the Euclidean space-time on which NC gauge theory 
is realized. 

In this paper, we investigate the vacuum expectation value (one point function) of these 
operators in detail. We develop an efficient procedure to calculate them through 1PI dia- 
grams. We perform explicit calculations on S 2 and S 2 x S 2 in various matrix models up to 
the two loop level. We identify the large N scaling behavior of the correlators to all orders by 
power counting arguments. We find that the correlators on S 2 and S 2 x S 2 exhibit different 
scaling behavior due to different dimensionality. In particular < TrA^A u > on S 2 x S 2 
indicates a 4d fractal like structure of it in contract to S 2 which remains to be semiclassical. 

The organization of this paper is as follows. In section 2, we investigate the effective 
action of quantum fields at the one loop level in preparation to calculate the correlators at 
the two loop level. We discuss how gauge invariance is respected in NC gauge theory. In 
section 3, we calculate the vacuum expectation value of this operator on S 2 up to the two 
loop level. We first evaluate it by the connected diagrams. We then evaluate it by 1PI 
diagrams after expanding A^ around the quantum solution which extremizes the effective 
action. We find that the both methods produce the identical answer as it can be argued on 
general grounds. The order parameter is found to consist of the Kronecker's 5 function of 
3d space in which S 2 extends. In section 4, we calculate the corresponding order parameter 
on S 2 x S 2 by 1PI diagrams. We find that it consists of the two independent tensors. We 
interpret the new tensor as the contribution from 4d fractals. In section 5, we explain that 
these operators also serve as the Wilson lines on homogeneous spaces [25] [26]. We compute 
the connected two point functions of these Wilson lines. We conclude in section 6 with 
discussions. 
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2 Effective action 



In this section, we investigate the one loop self-energy in NC gauge theory on S* 2 . This 
investigation is useful to compute the gauge invariant correlators up to the two loop level. 
Since NC gauge theory is realized by matrix models with finite N, the gauge invariance is 
exactly maintained. Our first concern in this section is how the gauge invariance is reflected 
by the one loop self-energy of gauge fields. 
Let us consider IIB matrix model: 

Sub = -j 2 Tr{h [ A fl ,A u ][A^A"] + l -^[A^]). (2.1) 

We expand A^ around a background which represents a fuzzy S 2 : 

A^ = f{p^ + a^). (2.2) 

where a M denotes the gauge field on it. p^ can be identified with the angular momentum 
operators in the spin / representation where N = n(2l + 1) in the case of n identical spheres. 
/ is the scale factor of this background. 

At the one loop level, the effective action consists of the gauge sector (gauge field and 
ghost) and the fermion sector contributions. The both contributions must be BRS invariant 
and the fermion sector contribution must be gauge invariant in particular. The fermionic 
contribution to the one loop effective action is 

-hrrlogil + r-SA^) 
= ~\tt[T ■ 5A^\ + \tt\Y ■ SA^T • 5A^\ + • • • , (2.3) 

where SA^X = [a^,X]. In this expression, the first and second terms represent the tadpole 
and self-energy of gauge field respectively. 

The gauge field self-energy (two point function) in the above expression consists of the 
planar and non-planar contributions: 

+ \ * " l "' r "Fo r "f i « |a " >n " {2A) 

where the symbol tr evaluates the trace over spinor indices only and 

PfjYj imi = \Pfj,,Yj imi \, 
QpYj 2Tn2 = lPn,Yj 2m2 \, 



We have also introduced the following average: 

<a"\X\a v > p = J2 *^23**«"23, 

31 ,j3,ni2,m3 

<a»\X\a u > np = J2 ^32*^23, 

j2,j3,m 2 ,m 3 

*a"23 = E Ujim^rYj^Yj^Yj^. (2.6) 

jimi 

For simplicity, we restrict our consideration to U(l) gauge theory here as it is straightforward 
to generalize it to U (n) gauge theory. 
Under the gauge transformation: 



a^P^A+^A], (2.7) 



the self-energy (2.4) changes as 



v 



v 



- < A\trT ■ P——T l/ -^—\a u > 
2 1 T-Q T-W 

-- < a^trT.-^— T • P— '— |A > 
2 1 M T-Q T-R 

~2 < ^ tr T^Q T ^ aU >p ~2 < A l* r f7^K >p 
+ \ < a v \tr^-^V v \K > p +i < a»\tr^-^T u \A > p 



= "8 E TrYl m2 (^Y j2m2 )[a\A] 



32,m 2 



+ 8 E ^(^^-)^m 2 K,A], (2-8) 

where we have retained only the linear terms in a. The non-planar part is gauge invariant 
by itself at the linearized level. 

Since the fermionic contribution to the effective action must be gauge invariant, it is 
somewhat surprising to find that the gauge field self-energy is not gauge invariant by itself. 
The resolution of this puzzle is that the gauge invariance is restored by the presence of the 
tadpole. The tadpole (one point function) is 

--TrF-SA- 1 



2 T-Q 

= 8 V TrY ] (9jLY- )a fl 
° 1 ' 1 nm2\r)2 nm.2) a 



Q 2 

- 8 E T r(%lY j2m2 )Yl m y. (2.9) 



J2, m 2 

'Qj. 

]2, m 2 
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The tadpole changes under gauge transformation as 

8 E TrYl m2 (^Y J2m2 )[a»,A} 

-8 £ Tr(^Y hm2 )Yl m2 [a^ A]. (2.10) 

Since (2.8) and (2.10) cancel each other, the fermionic contribution to the one loop effective 
action is explicitly shown to be gauge invariant at the linearized level. The lesson we draw 
here is that the transversality of the gauge field self-energy could be spoiled in NC gauge 
theory under the presence of the tadpole. 

We next investigate the gauge sector contribution. In this section, we deform the IIB 
matrix model action by adding a Myers term in such a way that S 2 is a classical solution. 

S IIB + l -ft, vp Tr[A,.A v ]A P - (2.11) 

Since IIB matrix model is a large N reduced model of lOd super Yang-Mills theory, we also 
evaluate the gauge sector contributions in the large N reduced model of D dimensional super 
Yang-Mills theory with the identical deformation. 
Firstly we evaluate the tadpole as: 

(D-2)TrSA-Q-^, (2.12) 

where D — 1 and —1 are the gauge field and ghost contributions respectively. It precisely 
cancels the fermionic contribution (2.9) when the physical degrees of freedom coincide (in 
super symmetric reduced models). 

The gauge field sector contribution to the planar part of the gauge field self-energy is 
evaluated in Appendix A as 

+(D - 2)(Q^R U + QyR^ - QyQ^ - R„R^ + (Q 2 + R 2 ) 8^)) \a v > p , (2.13) 

where 5^ is Kronecker's S in the 3 dimensional sub-space in which S 2 extends. 

Let us consider the D = 3 case which is relevant to a recent Monte-Carlo investigation [22]. 
Firstly the one loop self-energy is of the same order with the tree action if we fix f A N in 
the large limit. It is because (2.13) is 0(1/N) due to the 6j symbols in the interaction 
vertices. Secondly the self-energy is certainly positive for gauge fields which carry small 
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angular momentum. In the case of n identical branes with U(n) gauge group, there are 
zero-modes Oq. They represent the relative center of mass positions of the branes. For 
zero-modes, the gauge field self-energy can be estimated as 

~ ^(jog(N/n)-8)tra%aZ, (2.14) 

where tr is over the n dimensional subspace of zero-modes. We find that there arises positive 
mass term for zero- modes in the large N limit. Therefore the degeneracy is lifted at the one 
loop level. 4 

The fermionic contribution to the planar part of the gauge field self-energy (2.4) can be 
evaluated as 

-y < a»\ Q^p(Q»R» + Q»R» ~ Q ■ RS^ - 2V + V + ie lxvp P p ) \a v > p . (2.15) 

After combining (2.13) and (2.15), the planar part of the one loop self-energy of gauge field 
in a deformed IIB matrix model is found as 

\ < ^ g^( 4p2 V - 4 ^ + 8(<V - ^pP P ) ~ 16(^ - Ku)) W u > P ■ (2.16) 

It is invariant under the linearized gauge transformation <5a M = P M A just like the self 
energy of gauge field in ordinary gauge theory. However it contains mass and Chern-Simon 
type terms. After inspecting (2.16), we note that the one-loop self energy is positive for the 
components of the gauge field in the subspace in which S 2 extends. However it is negative 
for the components in the extra-dimensions as long as their angular momentum is small. 
They are always interpreted as scalars from gauge theory on S 2 point of view. In the case of 
n identical branes, (2.16) is consistent with the one loop effective action for zero-modes [18]. 

3 Correlators on S* 2 

In this section, we evaluate the expectation values of the following gauge invariant operator 
in NC gauge theory on S 2 : 

< ^TrA,A„ > . (3.1) 

Firstly we evaluate it through connected diagrams expanding around a classical solution. 
We then evaluate it through 1PI diagrams after expanding A^ around the quantum solution 
4 With finite N, it could become negative leading to an instability [22]. 
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which extremizes the effective action. We find that the both methods give the identical 
result as it can be argued to be the case. Based on the explicit calculations up to the two 
loop level, we identify the large N scaling behavior of this correlator. We further argue that 
this scaling behavior is valid to all orders. 
At the tree level, (3.1) is evaluated as 

^Trp,p v = J -l{l + 1) V ~ J — V (3.2) 

It is consistent with the fact that the space is classically a fuzzy sphere of the radius fl. In 
the classical limit this operator measures the shape and extension of a space realized in a 
matrix model. It is hence a natural operator to measure the geometry of the space at the 
quantum level as well. 
At the one loop level, 

P n 2 21 2i + 1 2n 2 

< Tra^ >= £ ] 1± T) ^ ~ —lo g (N/n)5^. (3.3) 

In a D dimensional bosonic reduced model, we also obtain the following tadpole contribution 
due to the presence of the one point function (2.12): 



f — < Tr Vil a v + a^pv >= -^—^NS^. (3.4) 



It is of the same order with the tree contribution (3.2) if we fix n 2 / f A N while (3.3) is 
suppressed by log(N)/N 2 in the large iV limit. Thus the quantum effects tend to shrink the 
radius of the sphere. This effect is caused by the logarithmic attractive potential between 
eigenvalues in bosonic models. In fact n 2 / f 4 N cannot be too large since otherwise (3.4) 
overwhelms (3.2) leading to unphysical results. 

It is an indication of an instability of the sphere when the coupling n 2 /f 4 N is strong 
enough. When D — 3, the stability of a fuzzy sphere solution can be estimated by the 
following one loop effective action with < A^>= ftp^: 

(\P 4 ~ \^f)Nl{l + 1) + N 2 log(3 + N 2 log(N/n). (3.5) 

The quantum vacuum is determine by minimizing it with respect to f3: 

1 r(/3 3 - P 2 f)N 3 + N 2 ^ = 0. (3.6) 



2n 2y ' ' J ' (3 
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In the weak coupling limit {(3 — > oo), this equation reproduces the classical solution 
/3 — f. The deviation from the classical solution can be estimated from (3.6) to the leading 
order of n 2 / f 4 N as 

f n 2 

The order parameter is estimated as 

N Trp ^ = T^v 1 - 4 7^ + ' ") 6 ^ (3 - 8) 

which is consistent with (3.4). We further note that the minimum of the effective action no 
longer exists when 

n 2 2 9 

" >! (3-9) 



f A N 3 3 

which indicates the instability of a fuzzy S 2 solution beyond this coupling constant. 5 In fact 
this argument has been found to be in good agreement with a Monte Carlo simulation [22]. 

Here we return to a supersymmetric model: the deformed IIB matrix model (2.11). At 
the two loop level, we obtain from the one loop effective action (2.16) 

P 



< Tra^a v > 

n 3 1/ ~\ 



pN (P 2 ) 2 Q 2 R 2 



(3.10) 



where 



<X> p = ^23**123, 
ji,J2,j3,mi,m2,m3 

^123 = TrYj 3m3 Y j2m2 Y jimi . (3.11) 

The symbol n.p. denotes the non-planar contributions. Each contribution in (3.10) is sup- 
pressed by 1/N 2 compared to the tree contribution when we fix n 2 / f 4 N in the large N limit. 
Although there are no tadpoles at the one loop level in supersymmetric models, they do 
arise at the two loop level. The total tadpole contribution to this order parameter at the 
two loop level is evaluated in Appendix B as 

f 2 n 3 ~ 172 1 

— < Tr Vil a v + a^p u >= <J^— < p2Q2R2 > P + n.p.. (3.12) 



'This critical coupling was first estimated by D. O'Connor. 
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The existence of the tadpole indicates the shift of the classical solution due to the quan- 
tum effects. Such a shift should also be calculable by minimizing the two loop effective 
action around = j3p^. We can read off the relevant effective action from our previous 
results[18][19] as 

- $ (32 + 4(f) - 12(1 -£)-*(!-/) + 24(1 - I) < ^ > p + n.p., 

(3.13) 

where the second line represents the two loop contribution. The second term proportional to 
{f 2 //3 2 ) denotes the contribution from the Myers vertices. The factor (1 — f / (3) represents 
the modification of the gauge field propagator from the minimal one. In this expression, we 
have retained only the linear contributions with respect to this factor. The third and fourth 
terms come from gauge and fermionic contributions and the last term comes from the Myers 
and cubic gauge vertices. 6 

The deviation from the classical solution can be estimated by minimizing (3.13) with 
respect to (3 to the leading order of n 2 / f 4 N as 

f n 5 344 

= >r + n - P - <3 ' 14) 

We thus find that 

N TrP ^ = (1 " ( /W } < PW P ) 12^ " V + ( } 

which agrees with (3.12). 

Let us consider the large N limit where we fix n 2 / f 4 N as the 't Hooft couplings. They 
are appropriate 't Hooft couplings in generic NC gauge theories on S 2 . It is because we 
obtain a factor n/N at each loop due to 6j symbols in the interaction vertices. The other 
factor n// 4 is generic in U(n) gauge theory. In such a large N limit, the correlators scale as 
follows to all orders 

^TrA^A, ~ nU^. (3.16) 

It is because the quantum corrections are small except the one loop tadpole contributions 
in bosonic models as we have seen in this section. 



6 Through this investigation, wc have realized that the coefficient 12 in the third term of (3.13) was 
erroneously evaluated to be 10 in [19]. We have corrected this error in [21]. 
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We can explain this scaling behavior by repeating the same argument with[19]. In this ar- 
gument, the large N limit of the amplitudes are estimated by their field theory counterparts. 
It is because we obtain 2d gauge theory in a semiclassical limit. In our case, the correlator 
is found to be logarithmically divergent at the one loop level in (3.3). It corresponds to the 
logarithmic divergence of the two point functions in field theory. Since 2d gauge theory is 
super renormalizable, higher loop corrections in 1PI diagrams are finite. 

As for the tadpole contributions, they can be estimated by the effective action. The 
effective action is quadratically divergent at the one loop level in bosonic models. It is the 
reason why the one loop tadpole contribution (3.4) is of the same order with the tree estimates 
of the correlators in the large N limit. The higher order contributions are suppressed by a 
power of 1/N 2 due to super renormalizability of 2d gauge theory. 

With U(l) gauge group in supersymmetric models, we may be able to consider another 
scaling limit where 1/ f 4 N 2 is adopted as the 't Hooft coupling. It is because another power 
of 1/N can be gained due to the cancellations between the planar and non-planar amplitudes. 
In this case, we find the following scaling behavior 

^TrA^-NS^, (3.17) 

where it is again dominated by the tree contribution. It is because the one loop corrections 
are 0(log(N)) and the higher loop corrections are 0(1) by the power counting arguments. 
We therefore conclude that the order parameter on S 2 is always given by a unique tensor: 
Kronecker's 5 in 3d in which S 2 extends. We note that the behavior of the order parameter 
remains the same with the classical limit. 

Another conclusion we can draw in this section is that we only need to calculate 1PI 
diagrams in our evaluation of the gauge invariant correlators. In this procedure we expand 

around the quantum solutions of the effective action. In fact we can give a generic 
argument to justify this procedure as follows. 

In order to calculate the correlators of TrA^A u , it is sufficient to calculate the free energy 
of the deformed action as 

e~ F = I dAd*be- s -^ TrA » A " 



J dAdipt 



SF 

< TrA^Ay > = (3.18) 

As it is well known, the free energy and the effective action are related by the Legendre 
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transform: 



w[j] = r[0] + j-0, 

' - 

where denotes a generic field. We can hence calculate F = W[0] by minimizing the effective 
action The derivations with respect to rf v are equivalent to the insertions of TrApA v 
operator into the connected or 1PI diagrams respectively. 

4 Correlators on S 2 x S* 2 

In this section, we compute the vacuum expectation value of the gauge invariant operators 
TrApAv on a 4 dimensional homogeneous space: fuzzy S 2 x S 2 . The order parameters on 
such a space can be evaluated in an analogous way as on S 2 in the preceding section. However 
they exhibit a different scaling behavior in the large N limit due to different dimensionality. 
We adopt the strategy to compute the 1PI diagrams by expanding A^ around a quantum 
solution which extremizes the effective action in this section. 
The first option is to consider the following deformation: 

SnB + ^Uv P Tr[A^A u }A p , (4.1) 

where f^ up /f are the structure constants of SU(2) x SU(2). A fuzzy S 2 x S 2 is a classical 
solution of this model. At the tree level, the order parameter is evaluated as 

f 2 f 2 o f 2 N - 

j^Trp,p u = J -l{l + 1) V ~ L^-S^ (4.2) 

where N = n(2l + l) 2 . We here assume that the both S 2 are of the identical size. The 
symbol 5^ u = f^pafvpa/^f 2 is Kronecker's 6 in the 6 dimensional sub-space in which S 2 x S 2 
extends. (4.2) is consistent with the fact that the eigenvalues of matrices are distributed as 
S 2 x S 2 in the 6d space with the radius fl at the classical level. 
At the one loop level, we obtain 

rp . _ n 2 (2j + l)(2p+l) 

N <lTa ^ > ~ f 2 Nj- p j(j + l)+p(p + l)^ 

Tl 

~ —2log(2)5^. (4.3) 

Since we identify n 2 / f 4 N as the 't Hooft coupling, it is of the same order with the tree 
contribution (4.2). 
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In order to evaluate the two loop contribution, we first evaluate the one loop self energy 
of gauge field. This evaluation is reported in appendix C. The order parameter < TrA^A u > 
itself is quadratically divergent by power counting on 4 dimensional spaces like S 2 x S 2 . In 
this case, we can safely ignore non-planar contributions since they are suppressed by a power 
of N. The leading part of the gauge field self-energy (C.5) is 

\ < a "\Q^( 4P %» - AP ^) K >p ■ (4-4) 
The two loop contribution to the correlator follows from this estimate as 

P 

— < Tra^a u > 

n 3 2 - 1 

:(4c^ - -V) < D2 ^ 2 D2 > P , (4-5) 



f 6 N x 3 P*Q*R? 



where 



<X> p = ]T *I 23 ^*123, 

^123 = TrYj. im . A Yj 2m2 Yj imi TrYp 3q3 Yp 2 q 2 Yp iqi . (4-6) 

We also need to determine < >= (3p^ by minimizing the two loop effective action. 
We can again read off the relevant effective action from our previous results [19] [21] as 

(V-V/)^ 

- f 4 (lQF! + 8(pFi - 12(1 - f -)F$ + 32(1 - L)(F* - Ff) + 24(1 - L)Ffj, 

(4.7) 

where the first and second line represent the tree and two loop contributions respectively. 
As in S 2 case: (3.13), the second term proportional to {f 2 /(3 2 ) denotes the contribution from 
the Myers vertices. The third and fourth terms come from gauge and fermionic contributions 
and the last term comes from the Myers and gauge cubic vertices. 

We have ignored the one loop contribution 0(log(N)) in comparison to the tree and the 
two loop contributions of O(N) in the large N limit where n 2 / f 4 N is fixed. It is a natural 
't Hooft coupling in 4d NC gauge theory as well since we obtain a factor of n/N due to the 
Qj symbols in the interaction vertices at each loop. Since 4d gauge theory is renormalizable, 
the over all degree of divergence remains the same to all orders. 
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The deviation from the classical solution can be estimated by minimizing (4.7) with 
respect to f3 to the leading order of n 2 / f 4 N as 

!-| = "7^(96^ + 4F|). (4.8) 

Since Ff and F% are 0(1) in the large A" limit, the two loop corrections are finite and small 
as long as the 't Hooft coupling is small. We have further checked that the two loop tadpole 
contributions are consistent with this result in appendix C. 

Alternatively we can evaluate the correlators on S 2 x S 2 in IIB matrix model itself. In 
fact it has been found that a fuzzy S 2 x S 2 is a consistent solution of IIB matrix model at 
the two loop level[19][21]. Furthermore the most symmetric case where the both S 2 possess 
the identical radii is favored by the effective action. It constitutes another evidence that 4 
dimensionality of space-time may be dynamically explained in IIB matrix model [27] [28] [29] . 
In this case, < >= fp^ is determined by minimizing the two loop effective action. 
Expanding around the quantum solution, the order parameter is given by the total of 
(4.2), (4.3) and (4.5): 

+ j 2 2log(2)5, u - -^(4^ - l^Ff 
= VN(^5^ + \log(2)5^-\ 3 ^(5, u -^ I/ )Fi), (4.9) 
where A 2 = An 2 / 'f A N . 

However we have found A 2 = \/^F% ~ 0.55 at the extremum of the two loop effective 
action. Therefore the two loop corrections in (4.9) are not small compared to the tree 
and the one loop corrections. Since a naive perturbation theory is not reliable, we adopt 
a resummation procedure to use the one loop exact propagator for gauge field. In this 
resummation, we insert the one loop gauge field self-energy into the gauge field propagator 
to all orders. 

Let us investigate the one loop gauge field self-energy in (4.4) which is of the form 
(PV - PpPvMP): 

( . ) ± v (2j 2 + l)(2p 2 + l)(2j 3 + l)(2p 3 + l) 

W[Jl,Pl) /S 2 ^ 2 , P3 02(J2 + 1)+P2(P2 + 1))0303 + 1)+P3(P3 + 1)) 

where we refer [30] for 6j symbols. In the large limit, we can adopt the Wigner approx- 
imation of 6j symbols to estimate the above amplitude[21]. In such an approximation, it 
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reduces to the following integral 



u = —-— C dYdZdMdN ' 



fH 2 T\ 2 h (Y + M)(Z + N) 

1 



x 



V2XY + 2YZ + 2ZX — X 2 — Y 2 — Z 2 — XYZ 



x - 1 (4.11) 

+ 2MN + 2NL - L 2 - M 2 - N 2 - LMN v 1 

where X = j 2 /l 2 , L = p\/l 2 . We can argue that this approximation is exact in the large N 

limit except at vanishingly small X and L. However such a region does not concern us since 

it does not contribute to the quadratically divergent order parameter. 

Integrating Z and M, we obtain 

u = A 2 f 4 dYdN 
Jo 

1 

X 



VX 2 + Y 2 + N 2 + 2XN + 2YN - 2XY - XYN 

x - 1 (4.12) 

VL 2 + N 2 + Y 2 + 2YL + 2YN - 2LN - YLN v ; 

For small X, L, the self-energy logarithmically diverges: 

^^fl(^) (4-13) 

as in Figure 1 where we plot u as a function of X = L. A full X, L dependence of u> is 
plotted in Figure 2. 

In terms of u>, the one loop exact propagator to the leading order in 1/P 2 expansion is 

1 1 tr + ?&^-.. (4-14) 



P 2 \ -\- LU P 4 1+u 

The order parameter is given in this resummation scheme as 

1 f 2 

- 6A ^+ f 2 N j- p3 (j + 1)+ p(p + 1)1+ u(j,p) 

1 v (2j + l)(2p+l) u(j,p) , 
^ Hi -1 ' 



Qf 2 N % j(j + 1) + p(p + 1) 1 + P) 



/AU- 



(4.15) 



In our approximation, it reduces to the following integrals 

^^xTITT^Ei^ (4 - 16) 
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1 - 



1 ' ' ' ' ' ' ' 1 

0.5 1 1.5 2 2.5 3 3.5 4 



Figure 1: Plot of u(X, L = X) against X: 

Line (a) is ^ (X . L 2 =X) and line (b) is 0.981n(4/X) + 0.74. 




Figure 2: Plot of uj(X, L) against X and L 
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We can numerically estimate the above integrals with A 2 ~ 0.55 as 

^(0.48^ + 0.91^). (4.17) 

Although this estimate is not exact, it is a universal conclusion that the gauge invariant 
correlator < -^TrA^A^ > scales as VN in 4 dimensional homogeneous spaces S 2 x S* 2 in the 
large N limit. It is because we can estimate the scaling behavior of the correlators to all orders 
by the power counting arguments. In 4d, gauge theory is renormalizable. Therefore the over 
all degree of the divergence of the correlators remains the same to all orders. The scaling 
behavior of the order parameter (y/~N) is the reflection of the fact that it is quadratically 
divergent. Hence this scaling behavior is valid to all orders with our identification of the 't 
Hooft couplings. 

The order parameter consists of the two independent tensors 5^ u and 5^ u . The first part 
can be interpreted as a distribution of S* 2 x S 2 extending in the 6 dimensional space just like 
our background. Since it is 0(R 2 ) in term of its radius R, we may identify R ~ iV 1//4 from 
(4.17). The non-commutativity scale (volume of a quantum) is set by the 't Hooft couplings. 

The second term proportional to 5^ u solely arises from quantum fluctuations. The iden- 
tical N dependence indicates that it also represents a 4 dimensional distribution. However 
5^ implies that it extends in 10 dimensions. We interpret it as the contributions from a 
fractal whose fractal dimension is 4. A concrete example is the branched polymers which 
have appeared in a previous study of IIB matrix model around commutative backgrounds [27] . 
Furthermore we have argued that 4d fuzzy homogeneous spaces are smoothly connected with 
the branched polymers[19]. The large N scaling behavior of the gauge invariant correlator 
(4.17) is consistent with our viewpoint. 

5 Connected two point functions 

In this section, we investigate the connected two point functions of the gauge invariant 
operators j^TrA^A u . In particular we are interested in the operators (/i ^ v) with no 
vacuum expectation values. Let us consider the connected two point functions of them 
on S* 2 . It is most convenient to expand the functions (or matrices) in terms of spherical 
harmonics Y\ m - Since p M can be identified with F ljti , this operator can be regarded as the 
gauge invariant completion of Tr (Yi^a u + Y^a^), namely a Wilson line on S* 2 which carries 
the minimum angular momentum[25] [26]. 
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At the tree level, the connected two point functions can be evaluated as 




(5.1) 



At the one loop level, there occurs the following contribution 



P P 
< —Tra^a v —Tra p a a > 



n 2 ^ 2j + l / 




(5.2) 



This contribution vanishes in comparison to the tree contribution (5.1) in the large N limit 
where either n 2 /pN or 1/ f 4 N 2 is fixed. We do find nonvanishing contributions to this 
correlator at the one loop level through the quantum corrections of the gauge field propagator. 
In the case of the 3D bosonic matrix model, the one loop self-energy of gauge field with small 
angular momentum scales as in (2.14). With the self-energy insertion into the propagator, 
the two point function behaves as 



where A 2 = n 2 / f A N. We thus find that the two point functions NG 2 on S 2 exhibits the 
logarithmic scaling violation at the one loop level. Since we can trace it to the logarithmically 
divergent one loop mass term, these logarithms may be resummed as 



In the case of supersymmetric models, the corresponding logarithmic scaling violations are 
absent since quantum corrections are finite to all orders in the large N limit. 

We next investigate the connected two point functions of these Wilson lines on S 2 x S 2 
in IIB matrix model. At the tree level, 




(5.3) 




(5.4) 



P P 




(5.5) 
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At the one loop level, there occurs the following contribution 

f 2 f 2 
< —Tra^—Trapae > 

n 2 _ (2j + l)(2p 



/ 4 ^ 2 ^0'0' + i)+p(p + i)) ; 

2 

j^logiN/ n) [5^5 ua + <5 M A P J ■ ( 5 - 6 ) 



This term is suppressed by 1/N in comparison to the tree term (5.5) in the large N limit 
where n 2 / f 4 N is fixed. Therefore the nontrivial contributions in the large N limit also come 
from the quantum corrections of the gauge field propagator on S 2 x S 2 . 
With the one loop self-energy estimation (4.13), we obtain 

+ Siur8vp)- (5-7) 



12n 

We thus find the logarithmic scaling violation of the two point functions G4 on S 2 x S 2 also. 
We can trace them to the wave function renormalization of gauge field. It is therefore likely 
that these logarithms may be summed by the renormalization group. Let us contemplate to 
renormalize the gauge field as 

= + Zj yP a^ p Yj CED Y p , 
Z 2 p = l-\ 2 log((N/n)(j 2 +p 2 )). (5.8) 

It is then natural to renormalize the two point functions of the Wilson lines as 

= Z- 2 G 4 . (5.9) 

By assuming Gf possesses a finite large N limit, we obtain the following renormalization 
group equation for bare G 4 

N-^G 4 = -A 2 G 4 . (5.10) 

This equation suggests that G 4 scales as (N/n)~ x2 in the large N limit with A 2 being fixed. 

It is clear that we need to renormalize NC gauge theory on homogeneous spaces to 
determine the scaling behavior of the Wilson lines. Admittedly our investigation in this 
section is exploratory and a more detailed investigation on this problem is necessary. We 
also need to study the effects of the non-planar contributions. 
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6 Conclusions and Discussions 



In this paper, we have investigated correlators of TrA fJi A u in matrix models of homogeneous 
spaces such as S 2 and S 2 x S 2 . It measures the geometry of Euclidean space-time on which 
NC gauge theory is realized. Since space-time and matter are unified in matrix models, the 
structure of space-time in general receives quantum corrections. This feature also suggests 
a close connection between NC gauge theory and quantum gravity. 

We have investigated these correlators in matrix models up to two loop level. We have 
developed a procedure to evaluate them through 1PI diagrams. We have determined the 
large N scaling behavior of the geometric order parameters to all orders. 

On 2 dimensional spaces S 2 , we have found that the expectation values of them are given 
by a unique tensor just like in the classical limit. This fact shows that the eigenvalues of 
matrices are always distributed as S 2 . On 4 dimensional spaces S 2 x S 2 , we have found 
that these order parameters consist of two independent tensors. Although the first tensor is 
due to a distribution of S 2 x S 2 , the second tensor indicates a distributions of a 4d fractal. 
We conclude that fuzzy S 2 x S 2 acquires a 4 dimensional fractal structure due to NC gauge 
theory on it in contrast to fuzzy S 2 . 

The operators with vanishing vacuum expectation values are the Wilson lines which carry 
the minimum momentum. We have investigated the two point functions of these Wilson 
lines on S 2 and S 2 x S 2 . We have found logarithmic scaling violations (or finite anomalous 
dimensions) of the Wilson lines. We note that it is a suppression effect for low multipoles 
which might be relevant to such a tendency observed in the cosmic microwave background. 
We believe that our procedure developed in this paper should be useful to study more generic 
Wilson lines. We hope that such investigations will further elucidate gravitational aspects 
of NC gauge theory. 

In the context of IIB matrix model, fuzzy S 2 x S 2 solutions may be interpreted as 
metastable D3-branes. We are most interested in the ultimate ground state into which these 
branes eventually decay. This is an analogue of tachyon condensation problem of unstable 
D-branes in IIB matrix model[31]. We also hope to make progress toward understanding 
this fundamental question. 
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Appendix A 

In this appendix, we evaluate the one loop gauge field self-energy in D dimensional bosonic 
reduced models with Myers term. We list the contributions from different diagrams sepa- 
rately. We only show planar contributions since non-planar contributions can be obtained 
from them as in (2.6) 

• From cubic gauge couplings: 

+(D - lXQpRv + Q V R„) - (4P 2 + Q 2 + R 2 )5^) \a v > p . 

(A.l) 



• From quartic gauge coupling: 

^ir < ° M| g^( (g2 + R2)5 ^ k >p ■ (A - 2) 

• From ghost: 

-\ < ^\q^{Q^ + QvR,)\a v > P ■ (A.3) 

• From Myers terms: 



-5<«'I^VK>p- (A.4) 



where 8^ = e^ pa e upa /2 is the projector into the 3 dimensional sub-space in which S 2 
extends. 

From Myers term and cubic gauge coupling: 



12 _L 

2 <a l Q*R 2 



< alX \7W51 ie ^P PP \ aV> p- ( A - 5 ) 



In total we obtain (2.13) in section 2. 

We further list the planar part of the ghost and fermion self-energy at the one loop level 
for completeness. 

• The ghost self energy: 

-<*!q4^>,- (A6) 
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The fermion self energy: 



1 

l Q~R 2 ' 



(A.7) 

Appendix B 

In this appendix, we evaluate the tadpole contributions to the following correlator on S 2 : 

P 

— < Trp^a v + a^v > • (B.l) 

The relevant cubic vertices in a deformed IIB matrix model are 

f 4 Tr( %[a u ,a p ] + -#Vk\V] + P fl b[a fl ,c]). (B.2) 

By attaching the one loop self-energy of gauge field (2.16) to the cubic vertices, we evaluate 
the pure gauge field contribution to (B.l) as 

n 3 ~ / 12 _4 1 104 1 \ 

/6JV ""V < Q 2 R 2 >p 3 < P 2 Q 2 R 2 >p 3 < P 2 Q 2 R 2 >p > + n ' P " { } 

We note that the contributions from the longitudinal part of the gauge field self-energy (the 
second terms in (3.10) and (B.3)) to the gauge invariant correlator cancel out each other. 
The remaining contributions are the followings. 

• The ghost contribution 

n 3 ~ 1 1 

fwH < Q^ >P + ^ (R4) 

• The fermion contribution 

,3 



n 



? / 64 1 64 1 x . 

<M ~ t < aw >p ~T < T^m >p ) + n - p - (R5) 



pN ""V 3 Q 2 R 2 p 3 P 2 Q 2 R 2 
• The quartic gauge coupling contribution 



pN Q 2 R 2 
In total we obtain (3.12) in section 3. 



n 3 ~ 9 
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Appendix C 



In this appendix, we evaluate the planar part of the one loop gauge field self-energy on 
S 2 x S 2 . We subsequently carry out a perturbative evaluation of the tadpole contributions 
to the following order parameter on S 2 x S 2 : 



P 

— < Trpyfln, + > 



ji,Pi,m.i,qi 

^a^23 — a ^imipiqi^ r ^jarna^j2rn2^jimiTrYp 3q3 Yp 2 g 2 Y piqi . 

jimipiqi 

The fermion sector contribution has the form as the following 



16 

y <a 



n 

^P 2 R 2 



-if, 



'Q P {Q-R) + R P Q 2 Q P R 2 + R P {Q-Ry 



fiup 



Q 2 



+ 



R 2 



-if, 



'Q P {Q-R) + R P Q 2 Q P R 2 + R P {Q-R) 



+ 



Q 2 R 2 
i-((Q 2 i? 2 - Q 2 R 2 )5 pu + (Q 2 R 2 - Q 2 R 2 )S 



Q 2 R 2 

Q„Qu + R^Ru - 5 i3 (Q 2 R 2 + Q 2 R 2 + 2(Q-R)(Q-R 



(C.l) 



The contribution to the gauge field self-energy from the gauge sector is 

+(D - 2){Q^R U + Q V R^ - Q V Q P - R V R P + (Q 2 + R 2 )8^)) \a v > p , (C.2) 
where in this context 



(C.3) 



(C.4) 



where Q P ,R P and Q P ,R P denote the components of the operators in the first and second 
three dimensional subspaces. 5ij is the Kronecker's 5 in the remaining dimensions. 

We restrict our attention in the case of D=10 for simplicity. Using (C.4) we obtain the 
1-loop gauge field self-energy on S 2 x S 2 : 

1 



n 

2 <a ^Q 2 ~R 2 ~ 



AP p P v - AP\ U - 8 5 pu + Sif^P" 
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-mil 



'Q P {Q-R) + R P Q 2 Q P R 2 + R P {Q-R) 



pup 



Q 2 



+ 



R 2 



, R . f ( Q P {Q-R) + R P Q 2 Q P R 2 + R P (Q ■ R) ' 



Q 2 



R 2 



16 



■((Q 2 R 2 - Q 2 R 2 )5, V + (Q 2 R 2 - Q 2 R 2 )~5 pu 



Q 2 R 2 

QuQv + Rp.Ru - 5 i3 (Q 2 R 2 + Q 2 R 2 + 2(Q • R)(Q ■ R 



| ^ > p 



(C.5) 



Now we can evaluate the order parameter (C.l) at the 2-loop level. There are con- 
tributions from cubic gauge, ghost, fermion, quartic gauge couplings which correspond to 
(B.3)-(B.6). In total we obtain 



n 



where 



Fl 
G p 



Q 2 R 2 + Q 2 R 2 + 2(Q-R)(Q-R) 



(C.6) 



2 

/_L_ 

\P 2 Q 4 R 2 



P 2 Q A R 4 



{P ■ R)Q 2 + (P • Q){Q -R) + (P- R)Q 2 + {P ■ Q){Q ■ R) 



(C.7) 



Finally we obtain 



< Trp^a u + a p p u > 



n 3 6. 



pu 



The result (C.8) is indeed consistent with (4.8). 



(4F| + 96F 4 P ) . 



(C.J 
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